Pulsars as the Sources of High Energy Cosmic Ray Positrons 
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Recent results from the PAMELA satellite indicate the presence of a large flux of positrons 
(relative to electrons) in the cosmic ray spectrum between approximately 10 and 100 GeV. As 
annihilating dark matter particles in many models are predicted to contribute to the cosmic ray 
positron spectrum in this energy range, a great deal of interest has resulted from this observation. 
Here, we consider pulsars (rapidly spinning, magnetized neutron stars) as an alternative source 
of this signal. After calculating the contribution to the cosmic ray positron and electron spectra 
from pulsars, we find that the spectrum observed by PAMELA could plausibly originate from such 
sources. In particular, a significant contribution is expected from the sum of all mature pulsars 
throughout the Milky Way, as well as from the most nearby mature pulsars (such as Geminga and 
B0656+14). The signal from nearby pulsars is expected to generate a small but significant dipole 
anisotropy in the cosmic ray electron spectrum, potentially providing a method by which the Fermi 
gamma-ray space telescope would be capable of discriminating between the pulsar and dark matter 
origins of the observed high energy positrons. 
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I. INTRODUCTION 



Dark matter particles annihilating in the Galactic Halo are predicted to generate a number of potentially observable 
products, including gamma-rays, electrons, positrons, protons and antiprotons. In contrast to gamma-rays, which 
travel along straight lines, charged particles move under the influence of the Galactic Magnetic Field, diffusing and 
losing energy, resulting in a diffuse spectrum at Earth. By studying the cosmic ray antimatter spectra, balloon or 
satellite-based experiments hope to identify signatures of dark matter. 

The PAMELA satellite, which began its three-year mission in June of 2006, is designed to measure the spectra of 
cosmic ray positrons up to 270 GeV and electrons up to 2 TeV, each with unprecedented precision [l[. Recent results 
show that the ratio of positrons to electrons plus positrons (the positron fraction) in the cosmic ray spectrum appears 
to stop decreasing and begins to climb quite rapidly between 10 GeV and 100 GeV 0- A similar trend was in fact also 
indicated by earlier experiments, including HEAT Q and AMS-01 although with lesser statistical significance. 

The behavior in the positron fraction observed by PAMELA is very different from that predicted for secondary 
positrons produced in the collisions of cosmic ray nuclides with the interstellar medium. Barring systematics (e.g. 
problems in the positron/proton separation at high rigidity), the data appear to indicate the existence of additional, 
primary sources of high energy positrons, such as dark matter annihilations taking place in the halo of the Milky 
Way 0, B S B B E3l ■ ^ should be noted, however, that while the observed spectral shape can be easily accommodated 
with annihilating dark matter, the normalization of this contribution to the PAMELA data requires a somewhat large 
annihilation rate. Such a large rate could, in principle, result from the dark matter possessing an annihilation cross 
section in excess of the valuepredicted for a simple s-wave thermal relic (av ~ 3 x 10~ 26 cm 3 /sec), for example due 
to the Sommerfeld effect 0, [T3, [H| or in the case that the WIMP is of non-thermal origin [l2j ■ Alternatively, a large 
annihilation rate could be generated as the result of significant inhomogeneities in the spatial distribution of dark 
matter, such as the presence of a large, nearby dark matter subhalo [l3j|. In either case, such scenarios are somewhat 
constrained by observations of gamma-rays [l4| , antiprotons 0, [l2j and synchrotron emission [l2l , [l5| . 

The challenges involved in explaining the PAMELA signal with annihilating dark matter lead us to consider a less 
exotic and purely astrophysical explanation for the observed positron flux. Energetic electron-positron pairs can be 
in fact produced in astrophysical sources, the leading candidate sites being pulsars - rapidly spinning, magnetized 
neutron stars, which emit pulsed electromagnetic radiation, as observed from Earth. 
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In this paper, we explore the possibility that the positron fraction reported by PAMELA may be generated by 
mature pulsars. Gamma-ray pulsars are predicted to produce energetic electron-positron pairs with a harder spectrum 
than that from secondary cosmic-ray induced origin, leading to the possibility that such sources may dominate the 
cosmic ray positron spectrum at high energies. We calculate the spectrum of such particles from known local pulsars 
(Geminga and B0656+14), and from the sum of all pulsars distributed throughout the Milky Way. As found in 
earlier studies , we find that both local pulsars and the sum of pulsars distributed throughout the Milky Way can 
contribute significantly to the observed spectrum. At 10 GeV, we estimate that on average only ~20% of the cosmic 
ray positrons originate from pulsars within 500 parsecs from the Solar System. If gamma-ray pulsars are formed at a 
rate of ~4 per century in the Milky Way, we find that the observed flux of ^10-20 GeV positrons could be plausibly 
generated in such objects. Similar conclusions were derived in [lj], [l8|. Above ~50 GeV, however, the positron 
spectrum is likely to be dominated by a single or small number of nearby pulsars. If the high energy electron-positron 
spectrum is dominated by a single nearby source, it opens the possibility of detecting a dipole anisotropy in their 
angular distribution (see also [H|). We find that such a feature could potentially be detected by the Fermi gamma-ray 
space telescope (formerly known as GLAST) [2(3], thus enabling a powerful test to discriminate between the pulsar 
and dark matter origins of the observed cosmic ray positron excess. 

The remainder of this article is structured as follows: In Sec. [Til wc review the known properties of pulsars and 
consider them as sources of high energy electron-positron pairs. In Sec. IIIIl we consider the nearby pulsars Geminga 
and B0656+14 and discuss their potential contributions to the cosmic ray positron spectrum. In Sec. lIVi we calculate 
the expected dipole anisotropy from nearby pulsars and compare this to the sensitivity of the Fermi gamma-ray space 
telescope. We summarize and draw our conclusions in Sec. [V] 



II. PULSARS AS SOURCES OF ELECTRON-POSITRON PAIRS 

In both models of polar gap [U, [22j and outer gap [23[ , electrons can be accelerated in different regions of the 
pulsar magnetosphere and induce an electromagnetic cascade through the emission of curvature radiation, which in 
turn results in production of photons which are above threshold for pair production in the strong pulsar magnetic 
field. This process results in lower energy electrons and positrons that can escape the magnetosphere either through 
the open field lines [25[ or after joining the pulsar wind [ID]. In this second case, the electrons and positrons lose 
part of their energy adiabatically because of the expansion of the wind. The energy spectrum injected by a single 
pulsar depends on the environmental parameters of the pulsar, but some attempts to calculate the average spectrum 
injected by a population of mature pulsars suggest that the spectrum may be relatively hard, having a slope of 
~1.5-1.6 [l8|. This spectrum, however, results from a complex interplay of individual pulsar spectra, of the spatial 
and age distributions of pulsars in the Galaxy, and on the assumption that the chief channel for pulsar spin down 
is magnetic dipole radiation. Due to the related uncertainties, variations from this injection spectra cannot be ruled 
out. Typically, one concentrates the attention on pulsars of age ^10 5 years because younger pulsars are likely to still 
be surrounded by their nebulae, which confine electrons and positrons and thus prevent them from being liberated 
into the interstellar medium until later times. 

Still, some energetics considerations can be done with simple analytical models; this will also help the understanding 
of arguments developed in the next Section. The rate of energy injection from a single pulsar in the form of pairs is 
limited by its spin-down power (the rate of energy loss corresponding to the slowing rate of rotation) . Assuming that 
this is simply due to the emission of magnetic dipole radiation, the maximum rate of energy injection can be written 
as (see e.g. |24|): 



V = ~^f- « lO 3 ^?, J?? P- 4 erg s-\ (1) 

where B\2 = i? s /10 12 G is the magnetic field at the surface of the star, i?io = i? s /10km is the radius of the star and P 
is the period of the star in seconds. The period P (gyration frequency £1) increases (decreases) with time as a result 
of the spin-down, according to 



(1 + t/T ) 

where tq ~ 3c 3 1 / (B 2 R^fl 2 ,) , I — (2/5)M s R 2 is the moment of inertia of the star with mass M s and Slo = 2tt/Pq is 
the initial spin frequency of the pulsar and Pq is the initial period. Numerically, this yields: 
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r = 7.4 x IOX2 2 (i^) R wPo y^rs. (3) 



It follows that the upper limit to the rate of energy deposit in the form of electron-positron pairs is 

2™° TO 



c = mn = l-m 2 — — - — „-• (4) 



TO 

In terms of the total energy injected in a time t after the pulsar birth, 

E tot (t) = \ml--^ = 6 x 10 43 P - 4 <£i2*5 7 1 , erg, (5) 

where £5 is the time in units of 10 5 years. Therefore, the total energy that a mature pulsar (t 3> To) has injected in 
the form of magnetic dipole radiation saturates to 

Etat « i/ng = 2.2 x 10 46 ( -4^) i? 10 ^o~ 2 ^g. (6) 



1.4M 



In the same assumption of a mature pulsar, we also have that Qq r* ^(t/To) 1 / 2 , where f2 is the gyration frequency 
measured today. For instance, for the Geminga pulsar (P = 230 ms) t f» 370,000 years and r R* 10 4 years (using 
£12 = 1.6 and Rio = 1.5), one has ^o ~ 166 s -1 (Po w 40 ms). For these values of the parameters, the total energy 
output of the pulsar is rather large, E to t ~ 10 49 erg, which could easily account for the high energy positron flux. It 
is worth stressing, however, that only a small fraction of this energy will eventually end up in the form of escaping 
electron-positron pairs, and thus this number should be treated as an absolute upper limit on the pair luminosity of 
a single pulsar. Qualitatively, the combined effect of a declining absolute luminosity [Eq. and of an increasing 
escape probability conspire in singling out typical ages of ^10 5 years for the pulsars expected to contribute maximally 
to the positron flux. 

To proceed in a more quantitative way towards the calculation of the overall spectrum from Galactic pulsars, one 
needs to adopt a model for the e + — e~ acceleration and escape probability from a single pulsar with a given magnetic 
field, period, etc. and then integrate over a Monte Carlo distribution of these typical parameters in a Galactic Pulsar 
population. The resulting injection spectrum we adopt follows from such a calculation in Ref. [l8j : 



— ^ « 8.6 x 10 38 iVi 00 (S e /GeV)- 16 exp (-E e /80 GeV) GeV" 1 s _1 , (7) 
dE e 

where iVioo is the rate of pulsar formation in units of pulsars per century. This expression corresponds to an average 
energy output in electron-positron pairs of approximately 6 x 10 46 erg per pulsar, i.e. to efficiency < 1% compared 
with the upper bound derived above. In the following, we inject this spectrum according to the spatial distribution 
of pulsars given in Rcfs. [l8l. [26|. 

Once electrons and positrons are produced, diffusion in the Galactic Magnetic Field regulates their motion. Unlike 
previous approaches to the problem, mostly based on a simple implementation of the leaky box model, we calculate 
the effects of propagation by solving the transport equation for electrons, including synchrotron and inverse Compton 
scattering losses: 



d dn e 
dtdK 



V 



K(E e )x7 



d 



B(E e ) 



dEe 



Q{E e ,x) 



(8) 



with a free escape boundary condition at 4 kpc above and below the Galactic Plane. Here dn e /dE e is the number 
density of electrons/positrons per unit energy, K(E e ) is the diffusion coefficient and B(E e ) is the rate of energy 
loss. We adopt K{E e ) = K a (l + E e /(3 GeV)) 5 with K Q = 3.4 x 10 28 cm 2 /s and S = 0.6, and B(E e ) = -bE 2 e with 
b = 10 _16 GeV _1 s _1 . Q corresponds to the source term described above. 

In Fig. [1] we show the spectrum of positrons and the positron fraction resulting from the sum of all pulsars 
throughout the Milky Way. In the upper panels, we show results for different rates of pulsar birth (one per 10, 25, 
or 100 years). The dashed line represents the baseline result neglecting the contribution from pulsars, including only 
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FIG. 1: The spectrum of cosmic ray positrons (left) and the positron fraction (right) resulting from the sum of all pulsars 
throughout the Milky Way. Also shown as a dashed line is the prediction for secondary positrons (and primary and secondary 
electrons in the right frames) as calculated in Ref. [27| . In the right frames, the measurements of HEAT [j| (light green and 
magenta) and measurements of PAMELA 0] (dark red) are also shown. We have used the injected spectrum reported in 
Eq. (J7J. In the lower frames, the upper (lower) dotted line represents the case in which the injection rate within 500 parsecs of 
the Solar System is doubled (neglected), providing an estimate the variance resulting from the small number of nearby pulsars 
contributing to the spectrum. 



the positrons produced as secondaries in the hadronic interaction of cosmic rays. In the right frames, the positron 
ratio is obtained considering, besides secondary leptons, also the primary electrons accounted as in [27|, to ease the 
comparison with previous literature. In the right frames, the measurements of HEAT Q (light green and magenta) 
and the measurements of PAMELA || (dark red) are also shown. 

In the lower frames of Fig. [T] we show the positron spectrum and the positron fraction for Aioo = 4 if the injection 
rate within 500 parsecs of the Solar System is doubled (upper dotted curve) or neglected (lower dotted curve). This 
provides an estimate of the relative importance of average nearby sources compared to the contribution from more 
distant pulsars. We will discuss this issue further in Sec. IIIII 

Interestingly, the best fit to the HEAT and PAMELA data appears to be obtained for A^ioo = 4, namely about one 
pulsar birth each ~ 25 years. It is worth noting that this number is only slightly higher than the typical estimates of 
the galactic core collapse supernovae rate, from which pulsars are formed. This rate has been estimated in a variety of 
ways, including from the scaling of rates in external galaxies, from the measured gamma-ray flux from galactic 26 Al, 
from historical observations of galactic supernovae, and from empirical upper limits from neutrino observatories (for 
a review, sec Ref. [28]). Also note that since the primary electron flux is determined from a fit to the absolute flux, 
which has uncertainties as large as ±50% around 10 GeV (see the cosmic ray review in Ref. [I^), the best-fit value 
of A^ioo extracted from the ratio is affected by at least an error as large. Additionally, in principle our numerical 
results could be modified if a different normalization for the diffusion coefficient were chosen; yet, the constraint on 
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the confinement time inferred from measurements of the boron-to-carbon ratio does not leave much freedom in this 
respect, at least in the energy region of interest here. The slope of the diffusion coefficient in energy is also expected 
not to be a critical parameter, since the positron "excess" has thus far been detected over a relatively narrow region in 
energy, where 6 should be virtually constant. Future observations will determine whether the actual excess positron 
spectrum extends to energies beyond ~100 GeV. 

The cutoff in the positron spectrum derived in our calculations is solely the result of the corresponding cutoff in the 
injection spectrum shown in Eq. iJTj). This cutoff is determined by the details of the development of the electromagnetic 
cascade in the pulsar magnctosphcrc and, even more importantly, by the distribution of periods, magnetic fields, and 
radii of mature pulsars. The exact value of the cutoff energy should, therefore, not be considered to be a robust 
prediction of the theory, although it represents a good estimate of the order of magnitude of the cutoff energy. For 
instance, in Ref. [25| it is argued that the typical energy of electrons and positrons in the cascade associated with the 
polar gap is E e ss P GeV which, for a period of 200 ms, would yield a maximum energy of the resulting 

pairs of ~ 400 GeV. On the other hand, the description of Ref. leads to different maximum energies (~80 GeV) and 
even different scalings of the cutoff with the parameters of the pulsar. It is therefore not very difficult to accommodate 
moderately higher energy positrons in the pulsar scenario, if future observations were to show that this is the case. 
Having listed possible caveats to the above results, what is remarkable is that, without any fine-tuning of the models, 
typical predictions for both the spectrum and the normalization can explain the positron fraction data very well. 

At energies lower than ~10 GeV, pulsars are not expected to contribute any appreciable flux because of the very 
hard spectrum, compared to the spectrum of secondary positrons produced in hadronic interactions of cosmic rays 
diffusing throughout the Galaxy. These secondary electrons approximately reproduce the steep spectrum of the parent 
nuclei and at low energies dominate the observed positron fraction. Since the spectrum of positrons from pulsars is 
important only at relatively high energies, we have neglected here the role of solar modulation. 



III. NEARBY PULSARS AS A SOURCE OF HIGH ENERGY COSMIC RAY ELECTRONS AND 

POSITRONS 

In this section, following earlier studies [HI [l9| , we re-explore the possibility that an individual or small number of 
nearby pulsars dominate the cosmic ray positron spectrum within the energy range studied by PAMELA. As argued 
in the previous section, in order to contribute significantly such a pulsar can be neither too young nor too old. The 
electrons and positrons from very young pulsars (i.e. much younger than 100,000 yrs or t 5 <C 1) are expected to 
be confined by the surrounding pulsar wind nebula. Furthermore, electrons/positrons in the PAMELA energy range 
require ~ 10 4 (D/100pc) 2 years to diffuse a distance, D. The electrons/positrons from very old pulsars, in contrast, 
will have diffused over a large volume, diluting their contribution to insignificant levels. There are at least two pulsars 
that appear to satisfy these constraints, Geminga which is estimated to be 157 parsecs from Earth and 370,000 years 
old, and B0656+14, which is estimated to be 290 parsecs from Earth and 110,000 years old. 

The spectrum of electron positron pairs at Earth is again calculated by solving the transport equation, but this 
time for a single source. Moreover we consider the case of a bursting source, namely one in which the duration of the 
emission is much shorter than the travel time from the source. For the pulsars discussed above this seems indeed to 
be appropriate. For the simple case of a power-law spectrum of electrons and positrons, Q{E e ) oc E~ a , injected a 
time t ago at a distance D from the Solar System, the spectrum observed at Earth is given by: 



dN, _ Q(E.) „ ^ f_D_ . , D/Dj „, 



where the diffusion distance scale is given by: 



D dil (E e ,t) « 2 x K{E e )t l (1 ^J^r^ . (10) 

y (1 - d)E e /E max 

As a result of energy losses, the maximum energy that electrons can reach the observer with is approximately given 
by B m ax ~ l/(6t) w 3 i^ 1 TcV. Since energy losses for electrons with energy ~10-100 GeV arc not crucial over a time 
period of a few hundred thousand years, Eq. (jT0|) approximately reduces to Ddn(E e ,t) ~ -J 4K(E e ) t, corresponding 
to ^300-500 v^5 parsecs over the energy range of 10 to 50 GeV. 

As a first example, we consider the nearby (D ~ 157 pc) Geminga pulsar, having an estimated age of 370, 000 
years [l6l.[l9|. In order to take into account the uncertainties in the injection spectrum, we carry out our calculations 
by assuming that electrons and positrons are injected with a spectral index a = 2 and 1.5, and exponentially cutoff 
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FIG. 2: The spectrum of positrons (left) and ratio of positrons to electrons plus positrons (right) from the pulsar Geminga, with 
the dashed lines as in Fig. 1. In the right frames, the measurements of HEAT [3] (light green and magenta) and measurements of 
PAMELA (dark red) are also shown. Here we have used an injected spectrum such that dN e /dE e oc E a exp( — E e /600 GeV), 
with a = 1.5 and 2.2. The solid lines correspond to an energy in pairs given by 3.5 x 10 47 erg, while the dotted lines require 
an output of 3 x 10 48 erg. 



above 600 GeV. As a default quantity, we consider a total energy of 3 x 10 erg injected as electron-positron pairs, 
which constitutes a few percent of the total spin down power of the pulsar. Our results for Geminga are shown in 
Fig. [21 In the left panel we plot the positron spectrum from Geminga for a =1.5 and 2. Again, the dashed line 
presents the spectrum from secondary positrons alone. The right panel shows the positron fraction for the two values 
of a and for two values of the total energy injected in pairs. The lower of these values in our default choice, 3 x 10 47 
erg (solid lines), while the higher, dotted lines represent the approximate energy required to generate the entire flux 
of excess positrons from Geminga alone (3.5 x 10 48 erg). We thus conclude that if Geminga were to dominate the 
observed positron fraction at high energies, it would have to transfer on the order of ^30% of its spin-down power 
into electron-positron pairs. Such a high efficiency to pairs appears unlikely. The (probably) subdominant role of 
Geminga is not particularly unexpected, given its relatively old age. 

We now turn our attention to the case of the pulsar B0656+14. B0656+14 is considerably younger than Geminga 
(approximately 110,000 years old), has a period today of P = 390 ms, and a current spin down luminosity that is 
approximately the same as Geminga. The spectrum of positrons and the positron fraction from B0656+14 arc shown 
in the left and right panels of Fig. [3j respectively. The lines are labeled as in Fig. [2] Because of the younger age 
of this pulsar, the flux of positrons expected from B0656+14 is somewhat higher than from Geminga, despite being 
somewhat more distant (D — 290 pc). In the case with a = 1.5, the predicted positron fraction can fit the PAMELA 
data if it injected ~ 8 x 10 47 erg in electron-positron pairs. This appears still large, but less extreme than in the 
Geminga case. 

In Sec. HH we found that nearby pulsars (D < 500 parsecs) are likely to dominate the pulsar contribution to the 
positron spectrum, especially energies above ~ 50 GeV. More distant pulsars, however, are still anticipated to play 
an important role in the lower energy range of the PAMELA positive excess. In Fig. 01 we show a combination of 
pulsar contributions to the high energy positron spectrum and the positron fraction. In particular, we include the 
contribution from all pulsars more distant than 500 parsecs (using a rate of 4 pulsars per century, as shown in the 
lower frames of Fig. [l} and the contributions from B0656+14 and Geminga (using 3 x 10 47 erg in electron-positron 
pairs from each and a spectral index of 1.5). As can be seen in the right frame, such a combination can provide a 
good fit to the preliminary measurements of PAMELA and accommodates for a rising positron fraction even beyond 
~ 100 GeV. 

A comment regarding this result is in order. In describing the contribution from all pulsars throughout the Milky 
Way, we have adopted a spectrum with an index of 1.6 and an exponential cutoff above 80 GeV. In contrast, we have 
used a higher cutoff (600 GeV) and slightly harder slope (1.5) for the nearby B0656+14 and Geminga pulsars. These 
results do not contradict each other, for the average injected spectrum from Ref. [l8T ] results from the sum over a 
realization with wide variability of pulsar properties, including injected energy and spectral index. In this respect, the 
results presented here are only demonstrative of the fact that a significant role of nearby pulsars, while not required 
to explain present data, is consistent with them, in which case they should dominate the high energy tail. If this is 
the case, interesting observational signatures are possible, one of which is discussed in the next section. In no case the 
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FIG. 3: As in Fig. [2] but from the nearby pulsar B0656+14. The solid lines correspond to an energy in pairs given by 3 x 10' 
erg, while the dotted lines require an output of 8 x 10 47 erg. 




E e+ (GeV) E e (GeV) 

FIG. 4: The positron spectrum and positron fraction from the sum of contributions from B0656+14, Geminga, and all pulsars 
farther than 500 parsecs from the Solar System. 

high-energy spectra presented here should be considered as a robust prediction, since they depend crucially on the 
detailed spectral properties of B0656+14, Geminga or and other nearby, mature pulsars that contribute significantly 
to the high energy positron spectrum. 



IV. DISTINGUISHING BETWEEN PULSAR AND DARK MATTER ORIGINS OF HIGH ENERGY 

COSMIC RAY POSITRONS 

The positron fraction reported by PAMELA taken alone is likely insufficient to distinguish between dark matter 
and pulsar origins of this signal. In this section we discuss an additional measurement which may help to resolve 
this issue. In particular, even after the diffusive propagation of electrons and positrons from pulsars is taken into 
account, at sufficiently high energies a small dipolc anisotropy should be present in the direction of the dominant 
nearby source(s). In a very general way, the anisotropy associated with diffusive propagation can be written as: 

5 _ /max - Jmin = 3K\V(dN e /dE e )\ 

+ 7 min c(dN e /dE e ) ' v ' 

where V{dN e /dE e ) is the gradient of the electron/positron density. The measurement of such an anisotropy in a 
statistical significant manner requires a large number of electron/positron events. For example, in order to detect an 
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FIG. 5: The dipole anisotropy in the electron+positron spectrum from a source 110,000 years old at a distance of 290 pc 
(B0656+14-like) and from a source 370,000 years old at a distance of 157 pc (Geminga-like). In each case, we have normalized 
the energy output to match the PAMELA data and have used a spectral shape of dN e /dE e oc E^ 1 ' 5 exp(— E e /600 GeV). Also 
shown as dashed lines is the sensitivity of the Fermi gamma-ray space telescope to such an anisotropy (after five years of 
observation). The Fermi sensitivity shown is for the spectrum integrated above a given energy. 

anisotropy at the 2cr level, one needs to fulfill the condition S > 2v / 2(Acv£obs)~ 1 ^ 2 , where N cv is the rate of events 
detected per unit time above a given threshold and f obs is the observation time. 

In addition to studying the gamma-ray sky, the Fermi gamma-ray telescope will also be able to measure a flux of 
electrons (and positrons, though without charge discrimination) at a rate of approximately 3 x 10 electrons per year 
above 10 GeV [301. This implies that Fermi should be able to detect (at the 2cr confidence level) a dipole anisotropy 
in the electron flux above 10 GeV if 5 > 0.05% in one year or S > 0.03% in 5 years 1 . 

In Fig. [5] we plot the level of anisotropy expected for a Geminga-like and a B0656+14-likc pulsar if they are 
responsible for the majority of the observed positron excess. The two dashed lines show the sensitivities of Fermi 
to anisotropy at 95% confidence level and at 5tr confidence level, after five years of observation (integrated above 
a given energy). We find that Fermi should be capable of identifying a single local source (or multiple sources in 
the same direction of the sky) if that source injected the bulk of its electrons/positrons within the last few hundred 
thousand years (the B0656+14-like and Geminga-like cases correspond to injection 110,000 and 370,000 years ago, 
respectively). If only a fraction of the high energy positrons observed by PAMELA originate from a given nearby 
pulsar, the corresponding solid lines shown in Fig. [5] should be multiplied (reduced) by this factor. Also note that 
B0656+14 and Geminga lie in similar directions in the sky, so they are expected to contribute the same overall dipole 
anisotropy. 

Alternatively, if dark matter annihilations throughout the Milky Way's halo are primarily responsible for the excess 
in the high energy cosmic ray positron spectrum, a small dipole anisotropy in the direction of the Galactic Center could 
also be generated. Fortunately, both B0656+14 and Geminga are in approximately the opposite direction, allowing 
for a potentially unambiguous discrimination between these possibilities. In the special and relatively unlikely case 
that a nearby dark matter subhalo in the direction of B0656+14/Gcminga is responsible for the observed flux, it 
would be difficult to distinguish between pulsar and dark matter origins using this technique. 

If anisotropy studies should prove inconclusive in resolving this issue, other information could be inferred from the 
shape of the positron fraction and of the overall electron/positron spectrum. Peculiar shapes can result from the 
superposition of the overall pulsar spectrum plus local contributions (see, for example, Fig. 4 or Ref. [HI). Future 
studies of the electron and positron spectra at higher energies will be especially important, as the spectral cutoff in the 
pulsar case is typically expected to be smoother and less sudden than that predicted from annihilating dark matter. 
Furthermore, combining electron/positron measurements with those of antiprotons, antideutcrons and diffuse gamma- 
rays may prove useful in distinguishing between these possibilities. Population studies of pulsars in the gamma-ray 




1 Note that, extrapolating from the data reported in Q, PAMELA can collect ~ 10 5 useful electrons plus positrons events per year above 
1.5 GeV. For the range of interest here this is several orders of magnitude below Fermi and unlikely to lead to a meaningful constraint 
on the anisotropy of high energy charged leptons. 
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band by Fermi are also expected to refine theoretical predictions and shed light on this issue. 



V. CONCLUSIONS 



The results recently reported by PAMELA strongly indicate the existence of a primary source or sources of high 
energy cosmic ray positrons. This result is unexpected and very interesting, even if of purely astrophysical origin. 
Several papers have appeared recently which discuss this signal within the context of dark matter annihilations. In 
this article, we have instead explored the possibility that the observed flux of high energy positrons is the result of 
electron-positron pairs being produced in nearby and galactic pulsars. We find that pulsars throughout the Milky 
Way, and a small number of nearby mature pulsars, such as B0656+14 and Geminga, could each plausibly generate 
the observed flux of positrons. To normalize the overall flux, on the order of a few percent of the pulsars' spin down 
power is required to be transferred into the production of electron-positron pairs. The prediction in the case of the 
sum of all pulsars in the Galaxy appears somewhat more robust in that it relies on the average statistical properties of 
these astrophysical objects rather than on the specific characteristics of nearby pulsars. It is remarkable that in this 
case, reasonable values for the parameters can lead to a positron spectrum consistent with the observations. Also, 
a pulsar origin would naturally fit the absence of an excess in the anti-proton data [3l[, since differently from dark 
matter scenarios no hadronic cascades are associated with the production of pairs in the magnetospheres. 

If a single or small number of nearby pulsars dominate the high energy positron spectrum, an observable dipole 
anisotropy in the electron-positron spectrum may be present. In particular, we find that the Fermi gamma-ray space 
telescope would be able to detect the anisotropy generated by a few hundred thousand year old or younger source 
with greater than 2a significance. A ~100, 000 year old source, such as B0656+14, would produce an anisotropy that 
could be detected by Fermi with greater than 5a significance. This provides a valuable test for distinguishing between 
pulsar and dark matter origins of the observed cosmic ray positron spectrum, complementary to the spectral shape 
and additional information from antiprotons and gamma-rays. 
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